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APPROXIMATING ROOT OF AN EQUATION 
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Abstract. We generalize Ramanujan method of approximating the 
smallest root of an equation which is found in Ramanujan Note books, 
Part-I. We provide simple analytical proof to study convergence of 
this method. Moreover, we study iterative approach of this method 
on approximating a root with arbitrary order of convergence. 



1. Introduction. 

Ramanujan method of approximating the smaUest root zq of an equation 
of the form. 



(1.1) ^vlfez'^-l 



k=l 

is found in Chapter 2 of Ramanujan note books, part-I without proof. It is 
assumed that aU other roots of (1.1) have moduli strictly greater than \zo\- 
For z sufficiently small, write 

^ oo 

1 — V A, 7^ ^ 

It follows easily that Pi = 1 and 

n-l 

(1.3) P„-^AjP„_,, n>l 

i=i 

Ramanujan gives (1.3) and claims, with no hypothesis, that P„/P„+i ap- 
proaches a root of (1.1) [HIS]. In this present study, we give analytic proof of 
his method and generalize this to approximate a root of nonlinear equation 
with the arbitrary order of convergence. 
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2. Generalization of Ramanujan method. 

Suppose that / be a real-valued function, defined and continuous on a 
bounded closed interval [a, b] of the real line that satisfies f{a)f{b) < 0. 
Assume, further, that / is continuously differentiable on {a,b). We wish to 
find real number a G [a, b] such that /(a) — 0. Let z be an approximate of 
a that satisfies [2 — a| < 1 and f'{z) ^ 0. Define 

(2.1) P(z)/(z) = l, Ifz^a 

Now applying Leibniz's rule^ of n*'' derivative on (2.1) and after simplifi- 
cation, we find that 

(2.2) P„(^) = --^ 5] (^^jP,(z)/("-'=)(z), n > 1 

where Pfc(z) is fc*'' derivative of l//(z), (fc = 1, 2, . . . , n) and Po{z) = 1/ f{z). 
Using (2.2), the next three values of P(z)'s are listed as follow as 
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P(z)- /^'H^),2^'W' 

p., f^'Hz) , j'iz)f'^'Hz) r{zf 

and so on. Now define for n > 1 
(2.3) Hn{z)- ^^^'^^ 



The operator I?" in (2.3) denotes n*'* order differentiation. Assume that, a 
is a zero of order one for f{z) then /(a) = and /'(a) 7^ 0. Set 

(2.4) Q{z) = ^ 



It follows easily from Taylor's series [3] that Q{z) — for z 7^ a and 
Q{a) = 1//'(q!). Since a is zero of order one of /(z), Q(z) is differentiable 
at a. Subtituting (2.4) in (2.3), we find that 
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Expanding Q{z) by Taylor's series at z = a on both numerator and denom- 
inator and applying the operator 13" and after simplification, we have 



Thus, we obtain 

(2.5) i?„(z)=0((z-a)"+i) 
Since jz — a| < 1 and letting n oo, we obtain 

(2.6) limi7„(z)=0 

n-i-O 

On the other hand, using Leibnitz rule on (2.3), we find that 

Since -D"(l//(z)) = Pn{z) and after simplfiiction, we obtain 

(2.7) H„iz) = z-a + n^-!^^ 
Letting n — >■ cxd and using (2.7), we obtain 



(2.8) q; = z + lim n- 



n— >oo 



Pniz) 



Remark 2.1. If we replace Pn{z) by n\Pn, Pn~i{z) by {n — l)!P„_i and 
assuming z — and /(O) = —1 in (2.1), (2.2) and (2.8), we get Ramanujan 
method of approximating the smallest root of an equation of the form (1.1) 

m- 

The recursion formula given in (2.2) may lead numerical instability in 
machine computation due to divisibility of f{z) when f{z) — ^ 0. To avoid 
such instability, multiply f{z)" on both sides of (2.1), (2.2) and set Tk{z) — 
Pk{z)f{z)''+\ we obtain To{z) = 1 and 

ri-l , ^ 

(2.9) T„(z) = ~Y^^ (^^jr,(z)/(z)"-'=-V'"~'^(^), n > 1 
Also, the equation (2.8) becomes 



(2.10) a = z + lim n/(z)%^, n = 1, 2, . . . 

n^oo ln(Z) 
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Thus, n*'' convergent of (2.10) gives 

OiK. z^nf {£)——-, n=l,2,... 

3. Some rational approximations. 

Ramanujan gives six examples to illustrate performance of his method to 
approximate the smallest root of equations. The approximations are given 
in the rational form. In particular, he provides an interesting rational ap- 
proximation that log 2 = 375/541. In this section, we present such rational 
approximations to m*'* root of some integers and logarithm of some rational 
numbers. 

3.1. m*** root of a rational number. Let us consider f{z) = — a and 
c be an integer which is very near to m*'' root of a. Now, using (2.2), we 
obtain the following recursion 

(3.1) r„(c) = - 1 (^) 5^^^^n(c) (c- - «)"-'-' c--' 

and 

(3.2) ^ = + 71(0-- a) 

where r = n— \\in<m, and r = m if n > m. The sequence of rational 

approximations of -\/2 are 

99 239 577 1393 3363 8119 19601 47321 
70' 169' 408' 985 ' 2378' 5741' 13860' 33461 

which is corected up to 3, 4, 4, 5, 7, 7, 8, 9 decimal places respectively. Table 

1 lists rational approximations of some irrational numbers using MATLAB. 



Table 1 . Rational approximations of Irrational numbers. 



S.NO 


Irrational No. 


Ral ioual Approx. 


siguilicaut, digits 


1 

2 
3 
4 


v/3100 


50623 
24337 
4603 
2302 
6064 
4813 
3110 
623 


10-digits 

9-digits 
8-digits 
7-digits 



3.2. Some Logarithmic values. Let us consider f{z) = log(2-|- 1) — a and 
= be an initial approximation (where a is a rational number). Table 2 
lists rational approximations of logarithmic numbers using MATLAB. 
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Table 2. Rational approximations of logarithmic values. 



S.No 


Logarithm Nos. 


Rational Approx. 


significance digits 


1 

2 
3 

4 


loge 1-5 
log, 2.0 
log, 3.0 
loge 1-2 


3858 
9515 
32781 
47293 
12667 
11530 
724 
3971 


7-digits 

6- digits 

7- digits 
6-digits 



3.3. Some other examples. Ramanujan lists the first ten convergents 
to the real root x + = 1, with the last convergent being 13/19 = 
0.684210 .... By Newton's method, the root is 0.682327804 . . .. Here, we 
provide three examples to approximate root of a;^ — 2a; — 5 = 0, — 3 = 
and X = sinx + ^ by taking the initial approximations xq = 2, xq = 1 and 
xo = 1 respectively. Table 3 lists first ten convergents of roots all the three 
equations. The roots of first two equations converge to 14 decimal places 
at 9*'* and 8*'' convergent respectively whereas the root of x = sinx + i 
converges slowly to 1.49730038909589 at 23*'' convergent. Wc sec from these 
examples, this method of approximating a real root is very slow process. 
However, we get good approximation when n is large. 



Table 3. Examples for Generalized Ramanujan method 



It 


x'^ -2x- 5=0 


- 3 = 


X = sin X + 





2.00000000000000 


1.00000000000000 


1.00000000000000 


1 


2.10000000000000 


1.10363832351433 


1.58288042035629 


2 


2.09433962264151 


1.09853245432531 


1.51838510578857 


3 


2.09455842997324 


1.09861223692174 


1.50077867834371 


4 


2.09455128205128 


1.09861230157476 


1.49783013943789 


5 


2.09455148653822 


1.09861228868606 


1.49735888023541 


6 


2.09455148143875 


1.09861228866513 


1.49730334991792 


7 


2.09455148154375 


1.09861228866810 


1.49729959647640 


8 


2.09455148154234 


1.09861228866811 


1.49730005778495 


9 


2.09455148154232 


1.09861228866811 


1.49730030987454 



4. Iterartive approach ON Generalized Ramanujan method. 

In this section, we study iterative approach on generalized Ramanujan's 
method to improve the accuracy of the root of nonlinear equations. If Zq be 
an initial approximate root of f(z) then using (2.10), 

(4.1) Zm+l = Zm + nf{zm) ^]!~^!^^'^ m = 0, 1, 2, . . . 
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where 

n-l 



(4.2) T„(z^) = -^ rjr,(z„0/(^m)""'"V^""'H^m), n>l 

k=0 ^ ^ 

Set Zfn^i — a =: Cm+i and z^^ ~ a — e^. Then by using (2.5) and (2.7), we 
find that the order of convergence of (4.1) is rt + 1. (i.e) 

(4.3) e„+i = O {e^+') 
Also, the condition of convergence is 

u A\ ,1 T„-i{z)Tn+iiz) 

(4.4) ?i + 1 — n „ , < 1 



T„(z)2 

Setting ri = 1 in (4.1) and after simphfication, we obtain Newton-Raphson 
method 

Zm+l^Zm~—, r, TO = 0, 1, 2, . . . 

f [Zm) 

Setting 71 = 2 in (4.1) and after simplification, we obtain Halley's method 

f{Zm)/ f'(.Zm) - n 1 9 

Zva+l — Z„i ^ \riii \i-cfr \9 TO — U, i,/,... 

By varying the values of n, we can find more formulas with the higher 
order of accuracy. Thus, this new iterative method is more generalization 
of Newton-Raphson method [3] and Halley's method [4 . 

Example 4.1. The equation f{z) = z — cosz has exactly one root a = 
0.73908513321516 between and 1 and starting with zq = 0, generate the 
sequence zi, Z2, ■ ■ ■ by taking to = 1 : 5 and n = 1 : 4 in (4-5) are listed in 
the Table ^. 



Table 4. Computation of root of z — cosz = 



m 


n = 1 


n = 2 


n = 3 


n = 4 


1 


1.00000000000000 


0.66666666666667 


0.75000000000000 


0.73846153846154 


2 


0.75036386784024 


0.73903926244631 


0.73908513352403 


0.73908513321516 


3 


0.73911289091136 


0.73908513321515 


0.73908513321516 




4 


0.73908513338528 


0.73908513321516 






5 


0.73908513321516 









We observed that from Table 4, when to = 1 and n = 1 : 4 (i.e. first 
row) the sequence corresponds to generalized Ramaujan method studied in 
section 2, which converges very slowly towards the root. Similarly, when 
n — 1 and to = 1 : 5 (i.e first column) the sequence corresponds to Newton- 
Raphson method, which converges to root at 5*'' iteration, second column 
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correspond to Halley method, which converges to root at 4*'' iteration and 
so on. 

5. Conclusion. 

In conclusion we note that the generalization of Ramanujan method of 

approximating the real root of nonlinear equations has been developed in 
this article. Firstly, wc have generalized Ramanujan's method of approxi- 
mating the smallest root of equations, which is found in Ramanujan Note 
books - I, to any real root with simple analytic proof. Secondly, we have 
proved that iterative approch of this method has arbitrary order of con- 
vergence. Moreover, we have shown that Newton-Raphson and Halley's 
method are special cases of this generalized Ramanujan method. 
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